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Abstract 

Recently it has been shown that the evolution of open quantum systems may 
be "unraveled" into individual "trajectories," providing powerful numerical 
and conceptual tools. In this letter we use quantum trajectories to study 
mesoscopic systems and their classical limit. We show that in this limit, 
Quantum Jump (QJ) trajectories approach a diffusive limit very similar to 
the Quantum State Diffusion (QSD) unraveling. The latter follows classical 
trajectories in the classical limit. Hence, both unravelings show the rise of 
classical orbits. This is true for both regular and chaotic systems (which 
exhibit strange attractors). 
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I. INTRODUCTION 



Quantum mechanics is nonlocal. Classical mechanics is local. Though it is widely 
believed that quantum mechanics is the more fundamental theory, attempts to describe 
classical phenomena by quantum equations are fraught with difficulties. Not only do the 
calculations become extremely cumbersome, but they are conceptually more difficult. 

There is no particular consensus on what it means to cross the "quantum — > classical" 
border. Many criteria have been suggested: rapid decay of macroscopic superpositions, 
localization in phase space, approach to coherent states, decoherence, near-determinism, 
positivity of the Wigner distribution, and non- violation of the Bell inequalities |]T]-§]. In 
fact, it is rarely necessary to choose; as a rule, all of these are satisfied for macroscopic 
systems. Both modern experiments and the emerging field of nanotechnology, however, 
are increasingly challenging this divide. As we probe the mesoscopic region, where both 
quantum and classical effects are important, it behooves us to have a better idea of what 
"classical" means. This is not only of theoretical interest; a better understanding should 
make it possible to produce more efficient numerical models of mesoscopic systems. 

In this letter, we show how the use of quantum trajectories, both continuous and discon- 
tinuous, illuminates the intermediate scales where neither a purely quantum nor a purely 
classical description is practical. These models not only simplify computation (one of their 
major original motivations), but describe when and how quantum nonlocality disappears 
in the classical limit. We illustrate our results by plotting regular and chaotic quantum 
trajectories, and show the limit where they recover classical regular and chaotic orbits. 

The Schrodinger equation is the basic dynamical law of nonrelativistic physics. However, 
strictly speaking, it applies only to the entire Universe as the only truly closed system. All 
other systems are open. It is well known that the environment is crucial for the emergence 
of classical features in a quantum system. The best-known example is the case of a mea- 
surement , leading to sharp values of some physical quantity. More generically, environments 
induce decoherence, which is closely related to the rise of classical properties Con- 
sequently, we shall concentrate on open quantum systems in the Markovian (time local) 
limit. 

Markovian open quantum systems are usually described by a master equation: 

p = -i[H,p\ +Y.{LmpLl - ]^{LlLm,p}^ (1) 

where p is the density matrix for the system, H its Hamiltonian, and the linear operators Lm 
describe the effects of the environment. However, equation (|I]) is not entirely satisfactory 
for our purpose. It describes only mean values, computed over both quantum and classical 
probabilities. Hence, density matrices by themselves do not tell us which features can be 
described classically and which require a quantum description 0. 

Let us illustrates this distinction for the mean values over a density matrix p of the 
position operator q and its square q^. These values give little idea of the actual degree of lo- 
calization of the particle. A large value of the spread (Tr(g^p) — Tr{qp)'^)^ could correspond 
either to delocalized particles (resulting from, e.g., the spreading of a single wavepacket) or 
to localized particles whose position is classically uncertain (resulting from, e.g., Brownian 
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motion). The evolution of localized particles can be efficiently computed with classical or 
semi-classical models, while delocalized particles necessarily require a less efficient but more 
complete quantum description. 

In the latter case, generally, there is nothing which corresponds to classical phase space 
trajectories. This is one of the chief difficulties in characterizing quantum chaos in a way 
similar to that used in classical dynamics. 

By unraveling the evolution of the density operator, one obtains, as we shall see, classical 
mixtures (i.e., classical probabilities) of quantum pure states (i.e. quantum probabilities). 
This allows one to distinguish quantum from classical and, at the same time, provides a 
powerful tool for practical computations. 

II. QUANTUM STATE DIFFUSION AND QUANTUM JUMPS 

In such an unraveling, one describes the system in terms of a normalized pure state 
\ip{t)) which follows a stochastic "trajectory" in Hilbert space. By averaging the pure state 
projector over all possible trajectories with appropriate weights, one reproduces the 

density operator p = M{\iIj) {ipD- This is analogous classically to replacing the Fokker-Planck 
equation for probability densities with a stochastic Langevin equation for single trajectories. 

Unfortunately, unlike the case of classical Brownian motion, the unraveling of the master 
equation ([l|) is not unique. Thus, there is some ambiguity in how one separates classical 
and quantum uncertainties, related to the ambiguity in identifying density matrices with 
quantum ensembles. In this section we consider two well known unravelings of (|1|). 

In quantum state diffusion (QSD), the (Ito) stochastic evolution equation for the nor- 
malized state vector \ip(t)) reads: 

mt)) = -tH\^{t))dt - I Y.(L]L, - 2{L])^L, + \{L,)^\'')\m)dt 

^ j 

+ Y.{L,-{L,)^)m))di, (2) 
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where the "noises" dE,j are complex-valued Wiener processes of zero mean M{dC,j) = 
and correlations M{d^jd^k) = 0, M{dE*d^k) = Sjkdt. This equation describes a continuous 
non-differentiable evolution similar to the familiar diffusive paths of a classical Brownian 
particle, but in Hilbert space instead of real space. QSD is the only continuous unraveling 
which satisfies the same symmetry properties as the master equation itself 0. 

Our second example is the quantum jumps (QJ) unraveling, which is closely related to 
photon counting. However, it can be defined for any Lindblad master equation P,^. The 
stochastic increment for the wave-function is 

|#(t)) = -tH\ijit))dt-WiL]L^ - {L]L,)^)\m)dt 

^ j 

The discrete Poissonian noises dNj assume the values or 1. Most of the time dNj = and 
the evolution is continuous and differentiable. However, whenever dNj = 1 there is a "jump" 
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to the state Lj\i/j{t)) / y {L^jLj)t. The dNj processes have mean values M|^) (dNj) = {L^jLj)^dt 
and correlations dNjdt = and dNjdNk = 6jkdNj. This means, essentially, that jumps occur 
randomly with an average rate (L^jLj). 

Let us illustrate these two unravelings for a simple example: the damped harmonic 



oscillator at finite temperature. H = ua^a, Li = ^Jrv)a^ and L2 = y {n + 1)70 where n is 
the thermal equilibrium mean photon number, n = {a^a)p, and 7 is the inverse relaxation 
time. For QSD one can show that any initial states tends to a coherent state: {ipit)) 
where a\at) = at|at) and at = {{q)ip + i{p)^)/ ^J{2). Furthermore, the evolution of at is 
governed by a classical equation, 

dat = —iujatdt — —atdt + ^/fryd^t (4) 

Hence, for this example at least, the QSD equation fully describes how the environment local- 
izes the quantum state down to a minimum Gaussian wavepacket, and how this wavepacket 
follows a classical trajectory. 

In [[7|, p!0| , p!T| it was argued that this is quite general, and not peculiar to the harmonic os- 
cillator. In ||12| the QSD equation was applied to the Kaotic Anharmonic OScillator (KAOS) 
system; in this case localization takes place despite competition with the delocalizing non- 
linear Hamiltonian. In general, when the system remains small (in the sense that it does 
not explore much of the phase space in units of h), the QSD trajectory presents no definite 
structure. However, when the parameters are such that the system explores a larger portion 
of phase space, the quantum trajectories exhibit a clear structure which approaches the 
classical strange attractor. 

Similar behavior has been demonstrated for other chaotic systems, including the weak 



link capacitor and quantum kicked rotor H], and the forced damped Duffing oscillator [13 



Once again, as we approach the classical limit, the structure of the strange attractor begins 
to appear. 

The case of QJ is quite different. No result analogous to that for QSD holds in the 
harmonic oscillator case. However, in |^ the same KAOS system was studied using the 
QJ equation, and again the classical strange attractor was observed whenever the system 
explored enough of phase space. 

This result was very puzzling at the time. We can now show that this behavior is also 
generic for QJ. In particular, we shall see that whenever the system is far from the origin of 
phase space (i.e., the harmonic oscillator ground state) its dynamics become very similar to 
QSD. 



III. LOCALIZATION OF QJ 

Consider a system in a state \ip{t)) with a Hamiltonian H and environment operators 
-^1 = s-iid L2 = ^/72a^. Define 

\m=Di-at)\m), (5) 

where at = (a)v(t) cind D{a) = expjaa''' — a*a} is the displacement operator. \4>(t)) is the 
state \ip{t)) displaced such that the mean position and momentum vanish: (0(t)|a|0(t)) = 0. 
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Note that \at\ measures the phase space "distance" of the state \'ip{t)) from the origin |0). 
Let Aa^ measure the width of the state ipt- 



Now consider the following limit. When the oscillator is well localized relative to its distance 
from the origin, so that \at\ » Aa^, then the rate of jumps {L'^L)t ~ \oit\^ is large, while the 
size of the jumps are small: ((ip{t)\{a^ — al){a — at)\ip{t)) = Aa^, and similarly for a"''. Hence, 
when the energy of the system is large relative to the "jumpy" evolution approaches 
closer and closer to a diffusion process like QSD, describing localized wavepackets following 
classical trajectories. 

Let us develop the QJ equation (|^) to first order in AQ;/|a|. To remove irrelevant phases, 
we use the 1-dimensional projector Pt = \'ip{t)){'ilj{t)\. From equation (|) one obtains: 



Using alV') = a\ip)+D{a)a\(t)) and aP^at = \at\'^Pt+al{a-at)Pt+CitPt{.a) -dD+O'^i^fj), 



{a}a)^ - {a})^{a)^ = {(f)(t)\a^ a\(f)(t)) . 



(6) 




(7) 



one deduces: 





and similarly for a^Pta. (These are our only environment operators.) Inserting this into (|^ 
yields 




+ (L,- - {L,)^)Pt 




\{L,)\ 



dW, 



(9) 



where we have made the approximation 




(10) 



with the dWj standard real Wiener processes {dWidWj = 6ijdt). Hence, 




(11) 
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Note that this equation (^TJ) preserves the norm of \ip{t)) and recovers (|l|) in the mean: 
p(t) = M{\ip(t)){ip(t)\) at all times. It is almost identical to the QSD equation (||), except 
that the noise is a real Wiener process multiplied by a phase which depends on the phase 
space position of the oscillator. This is the only choice of phase such that the state vector 
remains normalized in the diffusive limit. The QSD equation was chosen to be independent 
of the choice of phase. The solutions of equations (0) and ( pTD are very similar, as we shall 
see. 

We can see this by going to the classical limit and plotting the trajectories for the QJ 
equation. We have done this for an extension of the thermal model presented in the previous 
section. Note how in this case the jumps occur very frequently as we move away from the 
origin of phase space. (See figure 1.) This also works for the chaotic cases already studied 
with QSD; in the QJ well, we can see the emergence of classical structure. (See 

figure 2.) 

A similar derivation can be done for any choice of environment operators L which are 
linear in a and (e.g., x and p). 

IV. CONCLUSION 

We have seen how the use of stochastic evolution equations for pure states provides the 
ability to draw useful pictures of elementary quantum phenomena. These same equations 
provide powerful tools for practical numerical computations; in particular, one can exploit 
the existence of localized solutions to reduce the numerical difficulty of solving the equations, 
using the technique of moving bases p!^|-[T7[1 to produce more efficient computer simulations. 
Quantum unravelings have long suggested the possibility of separating quantum and classi- 
cal uncertainties for open systems; however, the ambiguity in the choice of unraveling has 
prevented any conclusions from being drawn as to the exact meaning of this separation. 

We now see that, as one approaches the classical limit, it should be possible to make 
this separation in a similar way for different unravelings. This suggests that, in a sense, the 
details don't matter: there is a single classical limit towards which a broad class of different 
"quantum trajectory" techniques all tend. In treating physical systems in the mesoscopic 
regime, this should make it possible to determine unambiguously which characteristics may 
be given classical and which quantum descriptions. This, in turn, may contribute greatly to 
the numerical solution of practical problems. 

We should also stress the limitations of this approach. In the first place, it depends on 
the "system/environment split," which will always include a certain ambiguity. While the 
exact boundary between system and environment is not important for macroscopic systems, 
for a microscopic or mesoscopic system it can be crucial. Second, we have limited ourselves 
to the Markovian approximation. This is mainly because little is known about the more 
general case. We expect that similar considerations apply to non-Markovian systems. 

We would like to thank Lajos Diosi, Francesca Mota-Furtado, Ian Percival and Riidiger 
Schack for useful conversations. This research was funded in part by the UK EPSRC, the 
EU Human Capital and Mobility Programme, and the Swiss National Science Foundation. 
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Figure 1. The harmonic oscillator with = 1 at finite temperature, the equilibrium 
displaced by a constant force /3; curves are plotted for the values /3 = 1,4, 10 and average 
thermal excitation n = 0.2. Each jump is marked; as the oscillator is displaced further from 
the origin, the jumps become frequent and small in effect, illustrating the transition from a 
"jumpy" trajectory to a diffusion process. 



Figure 2. The Poincare section of the forced damped Duffing oscillator in the chaotic 
regime. The value f3 gives the scale of the system, with (3 ^ oo the classical limit; this 
system classically exhibits dissipative chaos, and has already been investigated in this limit 
using QSD (see |]13[). We keep ^ = 1 constant. As the scale increases for j3 = 1,4, 10, the 



structure of the strange attractor clearly emerges, showing the diffusive limit of quantum 
jumps. This is similar to the behavior observed for QSD. The classical result is included for 
comparison. 
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